Electromagnetic Surface Modes at Interfaces with Negative Refractive Index 

make a "Not-Quite-Perfect" Lens. 
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Interfaces between media with negative relative refractive index generically support propagating 
electromagnetic surface polariton modes with large wavenumber. The relation of these modes to 
a recent prediction by Pendry of "perfect (real) image formation" by a parallel slab of negative- 
refractive-index material is analyzed. The "perfect image" theory is found to be incomplete without 
inclusion of a large- wavenumber cutoff that derives from a necessary wavenumber-dependence of the 
constitutive relations, and which controls the resolution of the image. 

PACS numbers: 78.20.Ci, 42.30.Wb, 73.20.Mf, 42.25.Bs 



Long ago, Veselago El noted that if an object could 
be viewed though a transparent slab of thickness w of a 
notional material with negative relative refractive index 
ri = — 1, simple ray optics shows that there is a focused 
image at a distance 2w in front of the object, which is real 
if the distance between the object and the back surface 
of the slab is less than w. In general, as seen below, 
the condition n = — 1 can be satisfied only for light with 
a frequency belonging to a discrete set of one or more 
special frequencies u>* that characterize an interface that 
supports "negative refraction" . 

Ray optics is only valid at lengthscales large compared 
to the wavelength of light, but recently Pendry B re- 
ported that if the condition n = — 1 can be supple- 
mented with the condition of perfect impedance-matching 
between electromagnetic waves in the two media at the 
special frequency ui* , the formal solution of Maxwell's 
equations, with a local effective-medium approximation 
for the constitutive relations, predicts that a real image 
formed by light at that frequency is "perfect", in that 
it reproduces the features of the object at all length- 
scales, however small. Pendry describes this counterintu- 
itive prediction from Maxwell's equation as "superlens- 
ing" ; superficially, his solution appears mathematically 
correct, but it has been controversial, and various com- 
mentators |j, [I| |(J have looked for flaws in his reason- 
ing. 

Ruppin |Q has found that, unlike conventional refract- 
ing interfaces, interfaces with negative refractive index 
support surface electromagnetic modes ("surface polari- 
tons"). In this Letter, I point out that Pendry's "per- 
fect image" result, and its limitations, can be understood 
from a degeneration of these modes in the impedance- 
matched limit. In the local effective-medium approxima- 
tion, the surface polaritons become degenerate and dis- 
persionless at the "lensing frequency" u)*, with no upper 
limit to their surface wavenumber. It is this unphysical 
feature that "explains" the "perfect image" prediction, 
but this does not appear to have been previously explic- 
itly noticed, either by Pendry (though he hints at 
a connection to "well-defined surface plasmons"), or his 



critics, or in Ref. [7J. The dispersionless surface modes 
are clearly a pathology of the approximation which ne- 
glects any wavenumber dependence of the constitutive 
relations. In fact, as with all real matter, the microscopic 
nature of the lensing medium will provide an "ultravio- 
let" (large wavenumber) cutoff. This cutoff is relevant for 
optics only in the "superlensing" limit, when it controls 
the actual resolution of real images. 

Since the direction of light rays in the ray-optics limit 
is the direction of the group velocity of the light, and the 
component of the wavevector parallel to the surface is 
conserved during refraction, it is easy to see that a surface 
with negative relative refractive index in some frequency 
range is an interface between media which both support 
propagating long-wavelength electromagnetic waves at 
those frequencies, but with group velocities of opposite 
sign. The group velocity of electromagnetic waves must 
be positive in both the low-frequency and high-frequency 
limits, but can be negative at intermediate frequencies, 
as shown in Fig. (Jlj) ; in an isotropic medium, transverse 
modes with negative group velocity must become degen- 
erate with a finite-frequency longitudinal mode as k — * 0. 

The Poynting vector E x H of a propagating electro- 
magnetic wave is parallel to the group velocity. Referring 
to the orthogonal triad of three vectors (E, H, k) that 
characterize propagating electromagnetic waves in an 
isotropic medium, Veselago Q introduced the term "left- 
handed media" to describe media with negative group 
velocity, as opposed to conventional "right-handed" me- 
dia with positive group velocity. (This terminology seems 
potentially misleading, as the "handedness" it refers to 
derives from the representation of the magnetic field as 
an axial vector, not intrinsic chirality of the medium.) 

Note that the model spectrum shown in Fig.(|l|) will 
only exhibit dissipation at the frequencies loti and uiT2- 
Some of the commentators [| have suggested that 
Pendry's calculation must fail in practice because it 
omits absorption effects, which, because the constitu- 
tive relations are frequency-dependent, must be present 
at some frequencies because of the Kramers-Kronig re- 
lation. However, this does not require dissipation in the 
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FIG. 1: In this model spectrum, the photon couples to the 
transverse components of two polarization modes (one elec- 
tric, one magnetic) in a material. The propagating modes 
at low and high frequencies are of the usual so-called "right- 
handed" type, with a positive group velocity, but are "left- 
handed" in the frequency range wti < oj < lul2, where the 
group velocity is negative. A surface has a negative refraction 
index in a frequency range where it is an interface between 
bulk regions with opposite signs of the group velocity. 



frequency range of interest. In fact, if the "left-handed" 
medium can be treated as loss-free, it must be a periodic 
structure, so a microscopic lengthscale is provided by its 
lattice spacing. If the material is a perfectly-periodic 
"photonic crystal", and possesses only one propagating 
mode (of each polarization) in the frequency range of in- 
terest, it can be treated as non-dissipative to the extent 
that non-linearity is negligible. 

In an isotropic medium with frequency-dependent 
but local (wavenumber-independent) constitutive rela- 
tions, the spectrum of electromagnetic waves predicted 
by the effective Maxwell equations is given by c 2 k 2 — 
w 2 e(w)/i(w), where e and /t are the dielectric constant 
and relative permeability. If waves in the two media have 
opposite group velocities, there will be a particular fre- 
quency oj* at which the wavelengths of propagating waves 
in the two media coincide: \\{oo*) = A 2 (cj*) = 2-K/k*] 
the condition n = —1 for a flat interface to produce a 
focused image is only realized at this special frequency. 
The group velocity v at this frequency is given by the 
expansion in Sou = oj — oj*: 



u} 2 e(w)fi(u) = c 2 k* 2 + 2c 2 k*v- 1 5uj + 0{5ojf 



(1) 



Let the interface be the plane z = 0. Following Ruppin 
I , I look for an "S-polarized" interface mode 



B z {x,y,z) 
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z < 0, 
z > 0, 



(2) 



where n\ and K2 are both positive; B z couples to H x and 
E v , and all are continuous at the interface: k\\E v = ojB z , 
and h\[ioH x = iaB z , where H 
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This only has a solution when \x\j\i2 is negative. Since 
by assumption ei/ii and 62/12 are both positive, 62/61 is 
also negative, and the interface has negative refractive in- 
dex. The source of these fields is an oscillating transverse 
surface polarization current 



J y (x,y) = J ( 



y i{k,,x-u}t) 
e ! 



k\\fioJo = ia(fii - ii 2 )Bq 



(4) 
(5) 



The condition giving the frequency of the mode comes 
from combining (0) with 
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C K 2 = C fc|| — U) 62/12- 

For frequencies close to lu* , the expansion 
for small positive 5k\\ = k\\ — k* , 5u>/6k\\ 
+ , where 

v _ V1V2 («i - «|) 
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can be used: 
v, as Skn — y 
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where V\ and t>2 are the group velocities in the two media 
at frequency oj* , which have opposite signs. 

This result is easy to understand: the group velocity 
of the surface mode is a weighted average of the bulk 
group velocities of the two media, with a larger weight 
given to the velocity in the medium with smaller k, into 
which the fields penetrate deeper. In the special case 
/ti = — /12, hi = K2, the competition between the two 
media is exactly balanced, and v vanishes; in this limit 
the predicted frequency of the mode becomes perfectly 
dispersionless with w(k^) = co* for all ku > k* . This 
corresponds to perfect impedance matching at frequency 
uj*: fii(uj*)/ei(uj*) — /i 2 (w*)/e 2 (w*), when no reflection 
of incident propagating waves with fcii < k* will occur. 

There is a second ( "P-polarized" ) mode, deriving from 
longitudinal surface polarization currents: 

J x (x,y) = Jg e il - k " x - ut \ (8) 
D z (x,y,z) = Dle^e^*-^, z < 0, 

= cg e -«2V(fc|| z>0, (9) 



where cj(/i2 — /i 
= ia'D z , with 
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= fc|| Jg , k\\UV = -luD z , and k^E* 
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Because ei(w*)/e2(w*) = /i2(^*)//ii(^*), when one of 
these interface modes has positive group velocity along 
the interface, the other group velocity is negative. The 
spectrum is schematically depicted in Fig. fj^) . 

It is now instructive to examine the coupling between 
the two sets of modes on the opposite faces of a slab of 
width d of a medium 2 embedded in medium 1 . The con- 
sistency condition (||) for the (B z , H x , E y ) slab polariton 



FIG. 2: (a): Generic spectrum of electromagnetic modes that 
propagate along a negative-refractive- index interface with sur- 
face wavenumber ku > k*: there are two surface modes, 
respectively with positive and negative group velocity; the 
shaded regions of the (uj,k\\) plane indicate where either one 
or both of the media supports propagating bulk modes, and 
uj* is the special frequency at which waves have the same 
wavelength X* = 2n/k* in both media, (b): The degener- 
ate spectrum predicted by the local Maxwell equations in the 
impedance-matched limit where "perfect lens" behavior has 
been predicted: the two surface modes become degenerate at 
the frequency uj* for all ku > k* . Such exactly-dispersionless 
degenerate surface modes are an artifact of the approximation 
which neglects wavenumber-dependence of the constitutive re- 
lations of the "left-handed" medium. 



modes becomes 



Hi = -^( tan h(K 2 d/2)) ±1 

Ml 



(11) 



where the ± distinguishes the modes which have sym- 
metric (+) and antisymmetric (— ) flux B z (z). When 
(fc|| — k*)d ^> 1, n 2 d is very large, and the splitting be- 
tween the symmetric and antisymmetric combinations of 
the modes on the two faces is very small, proportional 
to exp(— k\\d). In the other limit, they become the n 
— and n = 1 bands of modes confined to the slab, 
that emerge from the edge of the continuum of exterior 
propagating modes with initial group velocity v\ at fc|| ~ 
k* — n 2 ir 2 /2k*d 2 , n > (the n > 2 bands remain within 
the region of the (u>, ku )-plane where the slab medium has 
propagating modes) . Similar considerations apply to the 
(D z , E x . H y ) modes, but with fx replaced by e. The pre- 
dicted spectrum is shown schematically in Fig.(||). There 
is a single band-crossing of the surface modes (allowed 
because "S" and "P" polarizations do not mix) exactly 
at the frequency uj* , when K\ — k 2 — ko, at a surface 
wavenumber ko > k* given by: 



no- 



where 7 is the ratio 

> 2 (^*) + MjX^O 
^ 2 {uj*) - Mi( w *) 



7 



e - Kod = 7 2 < 1; 



djuj*) + e 2 (co*) 
ei(u)*) - e 2 (uj*) 



(12) 



(13) 



(the two expressions for 7 are equivalent because eifii 
= e 2 /i 2 at the frequency to*). Note that fco — y 00 i n the 
impedance-matched "perfect lens" limit where 7 = 0. 



FIG. 3: (a): Generic spectrum of the coupled electromag- 
netic surface modes of a flat slab of "left-handed" medium, 
with finite thickness d, embedded in a standard "right- 
handed" medium, calculated assuming local (wavenumber- 
independent) constitutive relations. The splitting between 
the even and odd combinations of the surface modes becomes 
exponentially small for ku large. An allowed band crossing 
occurs exactly at the frequency ui* , at fcii = ko (see text), 
(b): Predicted spectrum in the "perfect lens" limit (perfect 
impedance matching). The band-crossing point fco recedes to 
ku = 00, and for large fcii the surface mode frequencies differ 
from the "perfect lens" frequency uj* by exponentially small 
splittings proportional to exp( — fcyd). 



I now examine the solution of Maxwell's equations for 
the steady-state radiation field of an object illuminated 
with radiation at the special frequency uj* . The sources 
of the radiation field are the oscillating currents in the 
object that are excited by the illuminating field. I will 
assume the source current distribution is restricted to the 
region z < 0, and that the object is viewed though a slab 
of "left-handed" medium that is present in the region 
< Z\ < z < z 2 , where z 2 — z\ = d > 0. The source 
can be resolved into transverse Fourier components k|| 
in the x and y coordinates, and the radiation field of 
each Fourier component computed separately. Consider 



a Fourier component with kf. 



,*2 



k > 0, which 



produces an evanescent field in the region z > 0: 



B z = Bae^^-^'^Fiz,-/), 

k\\E v = uj*B z , -ik\\HoH x = /i 

D z = Doe^-^^Fiz,-^), 

k»H y = -uj*D z , -ih\e E x = e 



_j dB z 



(14) 
(15) 



where (using continuity of H x and E x at the interfaces) 

,-K\z-*A +/3e -«| Z - Z2 | ) ^g) 

' ' 1 (17) 

(18) 



F(z,7) = e~ KZ + ae: 



7 + 1 

7 2 — e~ 2 



When presented in this form, the solution of Maxwell's 
equations has a simple physical interpretation: F{z, 7) is 
the sum of the evanescent radiation field of the object, 
driven directly by the illumination, plus the radiation 
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fields of the induced surface polarization currents of the 
slab, driven by the radiation field of the object. The 
amplitudes of the two surface currents diverge when fey 
= fc ( ex P —2nd = 7 2 ), when the driving frequency lu* is 
exactly resonant with the coupled surface modes. This 
resonance recedes to fc|| = oo as the impedances are tuned 
to the "perfect lens" limit 7 = 0. 

When 7 = 0, the field of the small surface polarization 
current on the left interface at z\ exactly cancels the field 
of the object for z > z%, while for z < z\, it cancels 
the field of the second polarization current on the right 
interface at Z2- The field observed for z > z\ is then just 
the radiation field of the surface polarization current at 
Z2- If Z\ < d, this current has an amplitude that, as 
fcii — > oo, grows exponentially oc expfc||(d — z\) relative 
to the strength of the source radiation field. Since the 
object was taken to be to the left of the plane z = 0, 
the condition z\ < d is precisely the condition that a real 
image of some part of the object can be formed. 

The predicted exponentially-large amplification in the 
large-wavenumber limit limit occurs because the differ- 
ence between the coupled surface mode frequencies and 
the driving frequency lu* becomes exponentially small as 
fen — > oo (see Fig.(||b). It is a pathology of Pendry's so- 
lution that was noted in Ref. ^; however, while those 
authors recognized that, for \z — zi\ < d — z\, this ex- 
ponential amplification leads to an "ultra-violet" (large 
wavenumber) divergence of the expression for the pre- 
dicted radiation field of a point object at z — 0, they drew 
the incorrect conclusion that this divergence implied that 
the evanescent radiation field of such a point object can- 
not penetrate the slab. In fact, the divergence will be al- 
ways controlled by the short-distance cutoff provided by 
the physical nature of the interface; this may either be the 
maximum wavenumber of a surface polarization current 
{e.g., a surface Brillouin zone boundary), or a wavenum- 
ber at which fc-dependent constitutive relations move the 
surface mode frequencies away from near-resonance with 
the driving frequency lu* . 

Pendry's solution is only strictly valid without a cut- 
off in the case z\ > d, but remains non-singular in the 
marginal case Z\ = d, where the image of a source at 
z = is neither real nor virtual, but is exactly on the 
right surface z — Z2 of the slab. In this limit, there 
is an interesting interpretation of Pendry's formal re- 
sult: the replacement e — > — e, fi — ► —fj, is equivalent 
to the transformation (B, D, E, H) -> (-B,-D,E,H), 
which in turn is equivalent to time-reversal of the source- 
free Maxwell equations. The predicted "perfect image" 
formed by monochromatic radiation with frequency lu* 
is then analogous to the "spin-echo" after a 7r-pulse in 



magnetic resonance: the radiation field from a source at 
z = first propagates a distance d along the z-axis in 
the normal medium, with dispersion of its k|| -Fourier- 
components, then this dispersion is exactly reversed by 
subsequent propagation through an equal thickness d of 
the "time-reversing" medium to refocus the field at the 
"perfect image" point. However, this intriguing interpre- 
tation is spoiled by the inequivalence between a physical 
"left-handed" medium and an idealized "time-reversed 
vacuum" : the equivalence is only valid at a single fre- 
quency lu*, and as k — > the frequency of "left-handed 
light" approaches a finite value lul (see Fig.(ffl)); also, the 
medium will provide a large-wavenumber cutoff. 

In summary, I have shown Pendry's Q controversial 
"superlens" theory is incomplete, in that it fails to explic- 
itly include a large-momentum cutoff, which is needed to 
regularize the theory when it describes a real image. The 
high Fourier components of the image are produced by 
the near-fields of surface polariton modes @, || which in 
the "superlens" limit are exponentially close to resonance 
at the special lensing frequency u>*. The microscopic 
structure of the surface of the lensing medium will deter- 
mine this cutoff, which is unspecified in a local-effective 
medium approximation, and will limit the resolution of 
the image. While "not-quite-perfect" , it seems that this 
resolution can in principle be engineered to be signifi- 
cantly smaller than the wavelength of the illuminating 
radiation, without violating any fundamental physical 
principles. 

This work is supported in part by NSF MRSEC DMR- 
98094983 at the Princeton Center for Complex Materials. 
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